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Method and Apparatus for Sensing 
and Interpolating Color Image Data 

Field 

This invention relates to a method arid apparatus for sensing and 
interpolating image data, and more particularly, to an improved method and apparatus 
for efficiently and accurately sensing and interpolating color image data in an electronic 
imaging system. 

Background 

Electronic imaging cameras for recording still images are well known 
in the art. Such electronic imaging cameras may utilize two-dimensional electronic 
image sensing arrays such as charge-couple devices (CCD's) or complementary metal 
oxide semiconductor (CMOS) devices. Such two-dimensional image sensing devices 
typically comprise a regular array of image-sensing elements (cf. } pixels) that respond 
to incident illumination and provide an electronic signal corresponding to the intensity 
of this illumination. The electronic signals from the image sensing elements collectively 
form image data, which may be recorded in either analog or digital form. The signals 
may also, either immediately or subsequently, be displayed on viewing devices such as 
cathode-ray tubes or liquid crystal displays, or be printed to provide a hard copy. 

The image sensing elements, on their own accord, discriminate primarily 
the intensity of the incident light, and not its color. In order to record a color image, it 
is common to cover each element with a color filter chosen from a small collection of 
filter colors. 




Typical color imaging devices use one of two-color schemes. The first 
color scheme uses red, green and blue filters. The second color scheme uses cyan, 
magenta and yellow filters. In this way, each image-sensing element is made sensitive 
to the intensity of light having the color of its overlying filter. For example, an element 
5 covered by a red filter is responsive to red light, an element with a green filter to green 

light, and so on. 

The disposition of the color filters on top of the array of sensing 
elements can be one of several well-known patterns. One such disposition is a repeating 
pattern of red, green and blue vertical stripes, so that all sensing elements within a 

l 30 single column respond to the intensity of a single color of light. Other dispositions 

CO include checkerboard arrays. 

[j In such arrays, it is common practice to use a repeating pattern of colors. 

?i For example, one well-known and popular pattern of filters is known as the Bayer 

=== pattern. In the Bayer pattern, even numbered rows of the image-sensing device have 

ijl 5 alternating red and green filters, while odd numbered rows have alternating green and 

12 blue filters (the blue filters aligned with the green filters of the even numbered rows). 

\= This method of collecting color image data, although widely used, 

l± suffers from the problem that each image sensing element records the intensity of only 

a single color, (i.e., each image sensing element generates one-color image data) Image 
20 printing and display devices, on the other hand, commonly require information about 

the intensity of all of the colors at every pixel. The remedy to this problem is a process 
called "color recovery", in which the missing color intensities at each of the pixels are 
estimated from the data at surrounding pixels. 

Conventionally, the missing color data is filled in by 2-dimensional 
25 interpolation. For example, when a Bayer pattern is used, each blue-sensing element 

is surrounded by four green-sensing elements, and the missing value of green at the 
blue-sensing element may be estimated as the average of the data from the four green- 
sensing neighbors. 
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Conventional color recovery methods can produce images with 
objectionable artifacts, such as "color fringes", near sharp edges. The usual approach 
to solving this problem eliminates color fringes at the expense of image sharpness either 
by blurring the picture or by suppressing selected spatial frequencies, known as anti- 
aliasing. Treating the image in either manner has the disadvantage of producing a 
blurred image. 

A solution to the artifact problem is described by W.T. Freeman in U.S. 
Patent Nos. 4,663,655 and 4,774,565, both of which are herein incorporated. These 
patents set forth a color recovery scheme that uses nonlinear-filtering to produce an 
image that is sharp, correctly colored, and having a reduced incidence of the 
aforementioned artifacts. 

In the Freeman patents, a color-recovery process and an apparatus for the 
implementation thereof are disclosed. The process could be applied either on 1- 
dimensional data, produced by a single-line array of sensing elements with color filters, 
such as in a flat-bed image scanner, or on 2-dimensional data, produced by a 2- 
dimensional array of sensing elements with color filters, such as in an electronic 
camera. The method is easiest to describe in the case of a 1 - dimensional linear array 
having a repeating pattern of red-, green- and blue-sensing elements i.e., RGB RGB etc. 

First, in accordance with said method, for each color pixel type (e.g., 
red), there will be one actually measured value for that color at every third sensing 
element, with values for that color missing at the two intervening pixels. The first step 
thus is to fill in the missing red (R) values at green (G) and blue (B) elements by using 
any of a number of well-known interpolation methods. 

Secondly and similarly, the blue and green values are filled in for all 
pixels at which they have not been explicitly measured, using the same interpolation 
method. 

Third, at this point red (R), green (G) and blue (B) values exist for every 
pixel. One of these values is an explicitly measured value, and the other two have been 
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obtained by interpolation. The interpolated values are denoted by attaching a prime, as 
in R' or G 5 or B\ The interpolated data now have the form: 

/' 

Element 1 Element 2 Element 3 Element 4 Elements... 

R,G\B' R',G,B 5 R', G\B R, G', B' R\ G 5 B 5 

Fourth, at each pixel, the three color difference values, R-G, G-B, and 
f =i 5 B-R are calculated using either explicitly measured values or interpolated values such 

l=f as they are available at each pixel. 

L J Fifth, the R-G values across the row are then collected and subjected to 

a median filter of some predetermined length, typically an odd number in the range 5- 
^ 11- The same is done with the G-B values and the B-R values, so that at each pixel 

« 10 there is a complete set of median-filtered color-difference values. These are denoted 

IS as (R-G) mf , (G-B) mf and (B-R) mf . 

17- Sixth, at each pixel, the one explicitly measured color is combined with 

S3 the two median- filtered color difference signals that contain that color to produce 

improved estimates for the previously interpolated color values. For the example array 
15 described above, the first pixel has an explicitly measured value for R, but interpolated 

values for G' and B\ In that case, the R values are combined with the median filtered 
color differences (R-G) m f 3 and (B-R) mf to find the final estimates: 
G A = R - (R-G) mf 
B A = R + (B-R) m f 

where the A notation is used to denote a final recovered data value. This prescription 
is followed for each of the pixels, producing final color estimates that have the desirable 
20 properties claimed in the Freeman patents. 

The Freeman method may be practiced on two-dimensional color data 
as well, using essentially the same prescription. In that case, for each sensing element, 
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calculations are performed on a symmetrically disposed neighborhood of elements. In 
steps 1 and 2, standard 2-dimensional interpolation methods, such as bi-cubic 
interpolation, are used to produce the initial guesses for the missing colors at the 
element. 

5 Then, all of the color differences are found at each element, and in step 

5, a 2-dimensional median filter is applied to the neighborhood of the element to find 

the median-filtered color difference values. 

Finally, in step 6, improved estimates for the colors are calculated for 

colors that have not been explicitly measured at each pixel. This is done by combining 
C310 the one color that has been measured with all of the median-filtered color-difference 

values that contain that color, 
n It should be appreciated that in each case, any missing color is finally 

} J z determined by combining a measured color value with a single median-filtered color 

=2 difference. 

r=jl5 While the methods described in the Freeman patent provide desirable 

" results in certain applications, when practiced on comparatively large two-dimensional 

fy color data arrays, the number of calculations needed for complete color recovery 

becomes correspondingly large. Since the speed of the methodology is tied to the 

number of its required calculation, there is an existing need for a method for recovering 
20 missing color data in a two-dimensional color data array that provides results 

comparable to those obtainable through the Freeman methodologies, yet not requiring 

as many calculations. 

Summary 

25 In light of the above-mentioned need, the present invention provides a 

method (and apparatus) for recovering missing color data in a two-dimensional color 
array, said method (and apparatus) involving application of two one-dimensional non- 
linear interpolation processes. 
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In dealing with so-called real-time applications, it is sometimes found 
that the amount of computation involved in following known two-dimensional color 
recovery prescriptions becomes uneconomical to realize in special purpose circuitry, or 
too slow to implement as a computer program. By dividing the two-dimensional 
5 calculation into two one-dimensional calculations, fewer operations are required, yet 

good color fidelity is retained. 

In a preferred mode of operation, the central steps of the operation are: 
First, for each row (or other first dimension) of the image, apply a one-dimensional 
color recovery prescription to recover all of the colors present in that row, and then for 
Q 0 each column (or other second dimension) of the image apply another one-dimensional 

CO color recovery prescription to recover, at each element, any colors that were not already 

fj determined. The method, in greater detail, can be defined as electronically capturing 

H and processing image information by: (a) providing a two dimensional array of discrete 

image sensing elements, each discrete element capable of providing an electronic 
fgl5 information signal in response to incident illumination, said electronic information 

12 signal corresponding to the intensity of said incident illumination, each discrete element 

being specifically responsive to one of at least three predetermined colors; (b) 
I* obtaining first color image data by exposing the two dimensional array to image 

information bearing illumination such that each discrete element provides said 
20 electronic information signal, said first color image data comprising the collection of 

said electronic information signals; (c) recovering missing color information along a 
first dimension by (i) interpolating the first color image data along said first dimension 
to provide first interpolated color data for each of the discrete elements, (ii) forming a 
first difference channel between said first color image data and said first interpolated 
25 color data, and (iii) applying a first one-dimensional non-linear filter to said first 

difference channel, whereby a second color image data is obtained for each of the 
discrete elements; and (d) recovering missing color information along a second 
dimension by (i) interpolating the second color image data along said second dimension 
to provide second interpolated color data for each of the discrete elements, (ii) forming 
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a second difference channel between said second color image data and said second 
interpolated color data, and (iii) applying a second one-dimensional non-linear filter to 
said second difference channel, whereby a third color image data is obtained for each 
of the discrete elements. 
5 In light of the above, it is a principal object of the present invention to 

provide a method for recovering missing color data in a two-dimensional color data 
array. 

It is another object of the present invention to provide a method for 
recovering color data in a two-dimensional color data array, the method being 
"30 comparatively quick in its operation, yet providing good accuracy in respect of resultant 

53 color fidelity. 

LJ It is another object of the present invention to provide a method for 

/rj recovering color image data in a two-dimensional color data array, the method involving 

the conduct of two one-dimensional non-linear interpolations. 
£31 5 It is another object of the present invention to provide an apparatus well- 

u suited for conducting two one-dimensional non-linear interpolations for the purpose of 

^ recovering image color data in a two-dimensional color data array. 

M These and other objects of the present invention will become apparent 

from the following detailed description of presently preferred embodiments of the 
20 invention, taken in conjunction with the accompanying drawings. 

Brief Description of the Drawings 

Fig. 1 is a schematic block diagram of an apparatus for sensing and 
interpolating color image data (as in electronic imaging system) according to an 
25 embodiment of the present invention. 

Fig. 2 is a detailed schematic block diagram of the first color recovery 
filter 114 from Fig. 1. 
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Fig. 3 is a detailed schematic block diagram of a variant of the second 
color recovery filter 122 shown in Fig. 1, said variant being configured particularly for 
a so-called "3 -color" electronic imaging system. 

Fig. 4 is a detailed schematic block diagram of another variant of the 
second color recovery filter 122 shown in Fig. 1, said variant being configured 
particularly for a so-called "3-color" electronic imaging system. 

Fig. 5 is a detailed schematic block diagram of a variant of the second 
color recovery filter 122 shown in Fig. 1, said variant being configured particularly for 
a so-called "4-color" electronic imaging system. 

Fig. 6 is a more detailed block diagram of another variant of the second 
color recovery filter 122 shown in Fig. 1 5 said variant being configured particularly for 
a so-called "4-color 55 electronic imaging system. 

Figs. 7 and 8 are examples of configurations of discrete image sensing 
elements 20 employed in accordance with certain embodiments of the present invention. 

Detailed Description 

Conventionally, the recovery of non-sampled image data missing from 
the sampled one-color, or "first color", image data received from an electronic two- 
dimensional array of image sensing elements requires 2-dimensional processing in the 
region about each sensing element. The present invention can accomplish similar 
results — in respect particularly of certain array configurations (e.g., "Bayer"-type 
configurations) — with fewer calculations, by essentially "decoupling" the conventional 
2-dimensional process into two-incremental 1 -dimensional processes. The first 1- 
dimensional color recovery process generates intermediate second color image data 
from the sampled first color image data. The second 1 -dimensional color recovery 
process generates the final and desired third color image data from the second color 
image data. 

The variety of electronic image sensing arrays, for which application of 
the present invention is envisaged, is broad. Regardless, all will typically comprise a 
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two-dimensional array of discrete image sensing elements, wherein each discrete 
element is capable of providing an electronic information signal corresponding to the 
intensity of any illumination incident thereon, and wherein each discrete element is 
specifically responsive to one of at least three predetermined colors. 
5 As used herein, "color" shall mean a range of wavelengths within the 

electromagnetic spectrum that includes not only visible light, but the infrared and 
ultraviolet regions as well. The most common predetermined colors are either red, 
green, and blue, or cyan (a combination of the green and blue wavelengths), magenta 
(a combination of the red and blue wavelengths), and yellow (a combination of the red 
Q 0 and green wavelengths). The color sensitivity of the image sensing elements is typically 

CO accomplished by filtering input illumination so that different image sensing elements 

fj receive colored illumination. 

^ By exposing the two-dimensional array to image-information bearing 

« illumination, a collection of each electronic information signal received from each 

£□15 discrete element is obtained. This collection of signals forms the raw unprocessed one- 

12 color image data from which fully-recovered third color image data can be derived. 

The first step towards deriving the fully-recovered third color image data 
M is to first recover missing color information along a first dimension (e.g., along rows of 

the array). This is accomplished by interpolating the first color image data along the 
20 first dimension to provide a first-interpolated color data for each of the discrete 

elements, then forming a difference channel between the first color image data and the 
first-interpolated color data, and then applying a one-dimensional non-linear filter on 
the difference channel and combining with the first color data to obtain so-called first- 
recovered image data. The second color image data comprises a combination of the 
25 first-recovered image data and the first color image data. 

The second step derives fully-recovered third color image data from the 
two color image data by recovering missing color information along a second dimension 
(e.g., along columns of the array). More particularly, the second color image data is 
obtained by interpolating along the second dimension to provide second-interpolated 
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data for each of the discrete elements, then forming a difference channel between the 
second color image data and the second-interpolated data, and then applying a one- 
dimensional non-linear filter on the difference channel and combining with the second 
color data to obtain the so-called second-recovered image data. The third color image 
5 data comprises a combination of the second-recovered image data and the second color 

image data. 

As used herein, the term "difference channel" shall not be restricted to 
a channel containing a single color difference signal, but may include multiple color 
differences. For example, a single difference channel may contain both (red-blue) and 
(30 (red-green) color difference signals. Likewise, terms such as "first-interpolated data" 

f q and "second color data" shall be used to refer to data containing multiple color signals. 

" In general, an apparatus useful for implementing the inventive 

« methodology will comprise (a) a two-dimensional array of discrete image sensing 

elements, (b) a first color recovery module, and (c) a second color recovery module. 
fgl5 The two-dimensional array includes charge coupled devices (CCDs), 

" complementary-metal-oxide-semiconductor (CMOS) devices, and other known 2- 

ry dimensional sensor arrays in which individual sensing elements produce an output 

il signal that is responsive to the intensity of the illumination incident on the image 

sensing element. 

20 The discrete image sensing elements may be arranged in any 

predetermined pattern. The individual elements receive light reflected from an object 
and imaged through an optical system and, in response thereto, provides an electronic 
signal proportional to the intensity of the incident illumination. In one embodiment, the 
individual image sensing elements are arranged in a rectangular grid pattern, such that 

25 they form a plurality of rows and columns. Other arrangements are anticipated. For 

example, as shown in FIG. 8 and discussed in Example 3, arranging the image sensing 
elements according to a hexagonal grid pattern will result in a plurality of rows and 
columns in which the rows and columns are not perpendicular. (See, first and second 
dimensions 10 and 12 in FIG. 8). 
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The most preferred arrangement for the discrete image sensing elements 
involves use of alternating rows of red-green and green-blue filters that form a three- 
color checkerboard pattern known within the art as the "Bayer pattern". An example 
of a Bayer pattern is shown FIG. 7. As shown therein, one will appreciate that (a) no 
two contiguous discrete elements 20 that are along the first or second dimension 10, 12 
are specifically responsive to same color, and (b) no more than one discrete element is 
positioned contiguously between two discrete elements that are specifically responsive 
to the same color along the first and second dimensions 10, 12. 

As mentioned hereinabove, each discrete element is capable of providing 
an electronic information signal corresponding to the intensity of any illumination 
incident thereon, and wherein each discrete element is specifically responsive to one of 
at least three predetermined colors. When the array is exposed, the collection of signals 
received from each element is digitized and provided as first color image data. 

The first and second color recovery modules operate in predetermined 
separate dimensions of the first color image data to perform the two 1 -dimensional color 
recovery operations. For example, as explained in more detail below, if the individual 
image sensing elements are arranged in rows and columns, one of the color recovery 
modules may operate in a first dimension corresponding to the rows of data and the 
other color recovery module may operate in a second dimension corresponding to the 
columns of data. 

The first of the color recovery modules includes a first interpolator 
configured and arranged to interpolate in the first dimension the one-color image data 
received from the image-sensing array. The first interpolator provides output data that 
includes a plurality of first interpolated color data that corresponds to each individual 
image sensor element. Each of the first plurality of interpolated color data that 
corresponds to an individual image sensor element will be of a different sampled color 
than the first color data directly sampled by the individual image sensor element. 

The plurality of first color data and the plurality of first interpolated 
color data are combined and filtered in the first dimension to constitute a first-recovered 
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color data corresponding to each individual image sensor element. The output of the 
first color recovery operation is a plurality of second color data that includes the first 
color data and the first-recovered color data corresponding to each individual image 
sensor element. As used herein, first-recovered data refers to color data corresponding 
5 to each individual image sensor element that will be of a different color than that the 

first color data sampled directly by that particular image sensor element. For example, 
if an individual image sensor element provides first color data on red light incident 
thereon, the first-recovered color data could be green, blue or other predetermined 
colors other than the directly measured red color. 
30 The second of two 1 -dimensional operations is accomplished by the 

q second color recovery module, which includes a second interpolator configured and 

- t arranged to interpolate, along a different dimension, the plurality of second color data 

3 received from the first color recovery module. The second interpolator provides output 
z data that includes a plurality of second interpolated data corresponding to each 

4 5 individual image sensor element. 

The plurality of second interpolated data and the plurality of second 
LJ color data are combined and filtered in the second dimension to produce second- 

2 recovered color data associated with each image sensor element. As used herein, 

second-recovered data refers to a color data corresponding to each individual image 
20 sensor element that will be of different color than the first color data provided by the 

individual image sensor element and the first-recovered data described above. The 
output of the second color recovery operation is a plurality of third color data that 
includes the second color data and the second-recovered color data corresponding to 
each individual image sensor element. 
25 Although the present invention allows broad variation, to more tangibly 

illustrate its scope, particular embodiments and apparatuses for the implementation of 
the inventive methodology, and components thereof, are set forth in FIGS. 1 to 8. 

Figure 1 is a block diagram of color recovery apparatus 100 for 
recovering three colors of data. Color recovery apparatus 100 includes an electronic 
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imaging system (not shown) providing a plurality of first color data 1 12 as an input to 
a first color recovery module 1 14. The first color recovery module 114 operates on first 
color data 1 12 in the first dimension and provides an output of a plurality of second 
color data 128, comprised of the first color data and first-recovered color data. Second 
color recovery module 122 receives as an input the second color data 128 and provides 
as an output a plurality of third color data 132. 

First color data 1 12 is a schematic representation of the plurality of first 
color data. As described above, each of the first color data 1 12 corresponds to one of 
the individual image sensor elements in the image sensor array. 

In the first illustrative embodiment shown in Figure 1, data R 109 
represents a value of red filtered light measured at a particular individual image sensor 
element, data Gill represents a value of green filtered light measured at a particular 
individual image sensor element, and data B 1 15 represents a value of blue filtered light 
at a particular individual image sensor element. In the first embodiment, the first 
dimension of the first color data 112 corresponds to the rows of the plurality of 
individual image sensor elements. Similarly, the second dimension of the first color 
data 112 corresponds to the columns of the plurality of the individual image sensor 
elements. 

First color recovery module 114 includes first interpolator 116 and first 
color recovery filter 118. First interpolator 1 16 receives the plurality of first color data 
112 and interpolates this data in the first dimension to produce a first-interpolated color 
data 120. The (') notation indicates an interpolated color value. 

First- interpolated color data 120 and first color data 1 12 are filtered and 
combined in the first dimension by color recovery filter 1 1 8 to fully recover a second 
color data 128 for each individual image sensor element in the image sensor array, 
where the ( A ) notation indicates a recovered color data. 

Second color data 128 is provided to second interpolator 124 that 
interpolates the second color data 128 in the second dimension to provide second- 
interpolated color data 130. 
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Second-interpolated color data 130 and second color data 128 are filtered 
and combined in the second dimension by second color recovery filter 126 to fully 
recover a third color data 132 for each individual image sensor element in image sensor 
array. 

Figure 2 illustrates a more detailed description of the first recovery 
system 114 of the first embodiment of the color recovery system. In Fig. 2, Gl 
represents green color data contained in the rows of data containing red color data and 
G2 represents green color data contained in the rows of data containing blue color data. 
Two pairs of interpolators are used, one pair for the rows containing R and Gl data, 
and one pair for the rows containing the B and G2 data. The interpolator may be any 
conventional interpolator, i.e. linear, spline, sine, or other types as well. The individual 
choice of a particular type of interpolator is considered to be within the skill of the art. 

Interpolators 202 and 216 receive R data 218 and Gl data 220 
respectively and provide interpolated data via lines 205 and 215 respectively. In this 
embodiment, the R and Gl data alternate, and switch 204 is used to select and to couple 
the interpolated data on lines 205 and 215 to combiner 210 via line 209. Switch 206 
couples either R 218 or Gl 220 first color data to combiner 210. Switch 206 selects the 
other color from that of switch 204. In this embodiment, combiner 210 forms a 
difference data wherein the difference between R and Gl data corresponding to the 
same individual image sensor element is provided according to: 

(1) R-Gl',or 

(2) R'-Gl. 

This difference is provided via difference channel 21 1 to filter 212. 

Filter 212 is a non-linear filter that is provided to remove narrow peaks 
and valleys within the color-difference data. A particular class of non-linear filters 
capable of this type of response is rank-order filters, and in the preferred embodiment, 
a median rank-order filter is used. A suitable median filter is described in the U.S. 
Patent 4,802,108 assigned to Polaroid Corporation, which is herein incorporated by 
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reference. Combiner 214 combines the first color data selected by switch 206 with the 
filtered difference data 213. The combiner 2 1 4 provides: 

(3) R-F(R-Gl), or 

(4) G1+F(R-G1). 

FO represents the operation of filter 212. As can be seen from equations 
(3) and (4), the recovered green data G1 A are provided by equation (3) and recovered 
red data R A are provided by equation (4). These recovered color data are provided via 
line 217. 

Similarly, interpolators 222 and 236 receive B data 238 and G2 data 240 
respectively and provide interpolated data via lines 225 and 235 respectively. In this 
illustrated embodiment, the B and G2 data alternate, and switch 224 is used to select 
and to couple the interpolated data on lines 225 and 235 to combiner 230 via line 229. 

Switch 226 couples either B 238 or G2 240 first color data to combiner 
230. Switch 226 selects the other color from that of switch 224. The combiner 230 
forms a difference data wherein the difference between B and G2 data corresponding 
to the same individual image sensor element is provided according to: 

(5) B-G2*, or 

(6) B'-G2. 

This difference is provided via difference channel 231 to filter 232. 

Filter 232 can be a non-linear filter that is provided to remove narrow 
peaks and valleys within the color difference signal. As previously discussed, a median 
rank-order filter is used. Combiner 234 combines the first color data selected by switch 
226 with the filtered difference data 233. Also, combiner 234 provides: 

(7) B-F(B-G2), or 

(8) G2+F(B-G2). 

Where F() represents the operation of filter 232. As can be seen from 
equations (7) and (8), the recovered green data G2 A are provided by equation (7) and 
the recovered blue data B A are provided by equation (8). These recovered color data are 
provided via line 237. 
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The description of module 1 14 applies both to a 3 -color and to a 4-color 
system. In a 3-color system, Gl and G2 refer to the same physical color, whereas they 
will represent separate and distinct physical colors in a 4-color system. 

Figure 3 is a block diagram providing further detail for the second color 
recovery module 122 of the first embodiment of Fig. 1. This diagram applies 
specifically to recovery of colors within the odd numbered columns of the image for a 
3-color system. 

Interpolator 302 receives R data 301 that consists of first color data and 
interpolates this data in the second dimension to provide red data corresponding to 
every individual image sensor element. Similarly interpolator 314 receives B data 303 
that consists of first-recovered color data and interpolates this data in the second 
dimension to provide blue data corresponding to every individual image sensor element. 
As above, switch 306 selects interpolated or first-recovered blue data to provide as 
inputs to combiner 308, and switch 304 selects interpolated or first color red data to 
provide as inputs to combiners 308, 312 and 322. In this embodiment, first combiner 
308 forms the difference between the blue and red data corresponding to each individual 
image sensor element according to: 

(9) R-B\ or 

(10) R'-B A . 

This difference is provided via difference channel 309 to filter 310. 

Filter 310 can be a non-linear filter, such as those previously identified, 
that is provided to remove narrow peaks and valleys within the color difference signal 
309. Combiner 312 combines the red data selected by switch 304 with the filtered 
difference data 311. In this case, second combiner 312 provides: 

(11) R-F(R-B), or 

(12) R'-F(R-B). 

F() represents the operation of filter 310. Equations (11) and (12) 
illustrate that the second-recovered blue data B A are provided by equation (11). The 
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data provided by equation (12) are not used. These recovered color data are provided 
via line 3 17 to selector switch 328. 

A third combiner 322 combines the red data selected by switch 304 with 
green data 321. The green data may be either first color or first-recovered color data, 
5 depending on whether the image sensor element for which the computation is being 

made is an even numbered element or an odd numbered element in the second 
dimension of the image sensor. In either case, this data will be denoted simply as G. 
The third combiner 322 forms the difference between the red and green data 
corresponding to each individual image sensor element according to: 
Q0 (13) R-G, or 

£8 (14) R'-G. 

This difference is provided via difference channel 323 to filter 324. 
;i Filter 324 can be a non-linear filter that is provided to remove narrow 

« peaks and valleys within the color difference signal 323 and is also similar to those 

H5 previously identified. Fourth combiner 326 combines the green data 321 with the 

12 filtered difference data 325. In one embodiment the fourth combiner 326 provides: 

[U s (15) G+F(R-G),or 

JI (16) G+F(R-G). 

In this case, F() represents the operation of filter 324. As can be seen 
20 from equations (15) and (16), the second-recovered red data R A are provided by 

equation (16). The data provided by equation (15) are not used. These recovered color 
data are provided via line 327 to a third selector switch 328. 

As described above, this apparatus provides second-recovered blue data 
via line 3 1 7 for odd numbered elements in the second dimension of the image sensor, 
25 and it provides second-recovered red data via line 327 for even numbered elements in 

the second dimension of the image sensor. Selector switch 328, operating in synchrony 
with switches 304 and 306, transfers these recovered color data to line 329. 

The second color module 122 has been described above according to its 
operation for odd numbered columns along the second dimension of the image data. 
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The same apparatus can be used to perform the second color recovery for even 
numbered columns along the second dimension of the image data. The description of 
the apparatus is changed only in that the red input data R 301 and the blue input data 
B 303 are interchanged, and the three synchronous switches 304, 306 and 328 are 
initialized to positions opposite to those shown in Figure 3. In this configuration, the 
second-recovered color data on line 317 is red data, and the second-recovered color data 
on line 327 is blue data. 

The first and second color recovery modules illustrated in Fig. 2 and Fig. 
3 always generate recovered color values by using a combination of an explicitly 
measured color value with a filtered color difference. This prescription yields the best 
accuracy for the recovered color values. However, by relaxing this constraint and 
incurring a small penalty in color fidelity, we can substantially reduce the computational 
requirements of the second color recovery module. 

In this form of the color recovery method, which we shall refer to as the 
"abbreviated procedure", any missing colors recovered in the horizontal row operation 
are treated, in subsequent steps, as though they had been measured explicitly. The 
consequence of this change is that fewer filtered color-differences are needed in the 
vertical column operations, and therefore fewer of the computationally intensive 
filtering operations are needed. 

Fig. 4 illustrates another version of the second color recovery module 
122 in which the abbreviated procedure is implemented for color recovery in a three- 
color system. Interpolator 402 receives R data 401 that consists of both first and first- 
recovered color data and interpolates this data in the second dimension to provide 
second-interpolated red data corresponding to every individual image sensor element. 
Similarly, interpolator 414 receives B data 403 that consists of both first and first- 
recovered color data and interpolates this data in the second dimension to provide 
second-interpolated blue data corresponding to every individual image sensor element. 
Switch 404 selects different second-interpolated colors to provide as input to combiner 
408, and switch 406 selects different first color data or first-recovered color data to 
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provide as input to combiner 408 and combiner 412. In this case, combiner 408 forms 
the difference between the blue and red data corresponding to each individual image 
sensor element according to: 

(17) R-B\or 

(18) R'-B. 

The values of R and B (unprimed) may be either first or first-recovered 
color data. The difference is provided via difference channel 409 to filter 410. 

Filter 410 can be a non-linear filter that is provided to remove narrow 
peaks and valleys within the color difference signal 409 and it is preferred to be a 
median rank-order filter as previously discussed. Second combiner 412 combines the 
second color data selected by switch 406 with the filtered difference data 411. 
Combiner 412 in this embodiment provides: 

(19) B+F(R-B), or 

(20) R-F(R-B) 

where F(3 represents the operation of filter 410. As can be seen from equations (19) and 
(20), the second-recovered red data R A are provided by equation (19) and the second- 
recovered blue data B A are provided by equation (20). These recovered color data are 
provided via line 415. 

The above apparatus may also be modified to recover four colors per 
individual image sensor element. In the four color recovery system, the first color 
recovery module shown in Fig. 2 5 which recovers a second non-sampled color 
corresponding to each individual image sensor element, can be the same. However, the 
structures and apparatus of the second color recovery module must be more complex 
in order to recover the remaining two non-sampled colors corresponding to each 
individual image sensor element. 

The color data denoted as Gl and G2 are now understood to correspond 
to two different colors. For purposes of illustration, these colors will be described as 
two distinct shades of green, but they may be any two colors and this description should 
not be construed as limiting the scope of the invention in any way. 
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Figure 5 illustrates a second color recovery module designed for the odd 
columns in a four-color recovery system. In the embodiment illustrated in Fig. 5, Gl 
represents green color data contained in the rows of data containing red color data (odd 
rows) and G2 represents green color data contained in the rows of data containing blue 
color data (even rows). 

In the illustrated embodiment of Fig. 5, four (in two pairs) interpolators 
are used. The interpolation for each color is conducted in the second dimension, which 
in the illustrated embodiment is vertically along the columns. For each color either the 
first color data or the first-recovered color data is used in the interpolation. 

One interpolator will interpolate the B data, one interpolator will 
interpolate the R data, one interpolator will interpolate the G2 data, and one interpolator 
will interpolate the Gl data. Thus, all data associated with an individual image sensor 
element will now have four colors, a first color data, a first-recovered color data, and 
two second-interpolated color data. The interpolators may be any conventional 
interpolator known in the art. 

Interpolators 502 and 518 receive B data 514 and R data 516 
respectively, and provide second-interpolated data via lines 503 and 519 respectively. 

Switches 504 and 520 are used to select and to couple the data to 
combiner 506. The switches couple first color red data 516 and second-interpolated 
recovered blue data 503, or else first-recovered blue data 514 and interpolated red data 
519, based on whether an odd row or an even row of data is being combined. 

If an odd row is being combined, then R 5 1 6 data and interpolated values 
of B A 503 are used. If an even row is being combined B A 514 data and interpolated 
values of R 5 19 are used. Combiner 506 forms a difference data wherein the difference 
between R and B data corresponding to the same individual image sensor element is 
provided according to: 

(21) R-B A \or 

(22) R'-B A . 

This difference is provided via difference channel 507 to filter 508. 



20 



Filter 508 is a non-linear filter of the previous types disclosed. 
Combiner 510 combines the red data selected by switch 520 with the filtered difference 
data. In this case, combiner 510 provides: 

(23) R-F(R-B), or 

(24) R'-F(R-B). 

F0 represents the operation of filter 508. As can be seen from equations 
(23) and (24), the second-recovered blue data are provided by equation (23). The data 
provided by equation (24) will not be used. 

Similarly, interpolators 524 and 540 receive G2 data 536 and G1 A data 
538 respectively, and provide interpolated data via lines 525 and 541 respectively. 
Switches 526 and 542 are used to select and to couple the data to combiner 528. The 
switches couple interpolated G2 data 536 and first-recovered G1 A data 538, or else 
measured G2 data 536 and second-interpolated recovered G1 A data 541, based on 
whether an odd row or an even row of data is being combined. 

If an odd row is being combined, then second-interpolated values of G2 
525 and first-recovered values of G1 A 538 are used. If an even row is being combined 
then first color values of G2 536 and second-interpolated recovered values of G1 A 541 
are used. 

The combiner 528 forms a difference data wherein the difference 
between G2 and Gl data corresponding to the same individual image sensor element 
is provided according to: 

(25) G2'-G1 A , or 

(26) G2-G1 A \ 

This difference is provided via difference channel 529 to filter 530. 

As previously discussed, filter 530 can be a non-linear filter that is 
provided to remove narrow peaks and valleys within the color difference data. 
Combiner 532 combines the color data selected by switch 526 with the filtered 
difference data. In one embodiment combiner 532 provides: 
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(27) G2'-F(G2-G1), or 

(28) G2-F(G2-G1), 

where F() represents the operation of filter 530. As can be seen from equations (27) and 
(28), the second-recovered green data G1 A are provided by equation (28). The data 
5 provided by equation (27) are not used. 

A third combiner 544 receives red data selected by switch 520 and G2 
data selected by switch 526. These switches couple first color R data 516 and second- 
interpolated G2 data 525, or else first color G2 data 536 and second-interpolated R data 
5 1 9, based on whether an odd row or an even row of data is being combined. 
£310 If an odd row is being combined then measured values of R 516 and 

JO interpolated values of G2 525 are used. If an even row is being combined then 

n interpolated values of R 519 and measured values of G2 536 are used. 

H Combiner 544 forms a difference data wherein the difference between 

=5 R and Gl data corresponding to the same individual image sensor element is provided 

f45 according to: 

(29) R-G2', or 

ry (30) R'-G2. 

The difference is provided on difference channel 545 to filter 546. 

Filter 546 can be a non-linear filter that is provided to remove narrow 
20 peaks and valleys within the color difference data. The filtered difference data 

produced by filter 546 is provided on line 547. 

Combiner 548 combines the color data selected by switch 520 with the 
filtered difference data 547. In one embodiment combiner 548 provides: 
(31) R-F(R-G2) 5 or 
25 (32) R'-F(R-G2). 

In this case, F() represents the operation of filter 546. As can been seen 
from equations (31) and (32), the second-recovered green data G2 A are provided by 
equation (31). The data provide by equation (32) are not used. 
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Combiner 550 combines the color data selected by switch 526 with the 
filtered difference data 547. In one embodiment combiner 550 provides: 

(33) G2 , +F(R-G2), or 

(34) G2+F(R-G2). 

F0 represents the operation of filter 546 in these equations. The second- 
recovered red data R A are provided by equation (34) and the data provided by equation 
(33) are not used. 

The second color recovery module 122 has been described above 
according to its operation for odd numbered columns along the second dimension of the 
image data. The same apparatus can be used to perform the second color recovery for 
even numbered columns along the second dimension of the image data. The description 
of the apparatus is changed only in that: 1) the blue input data 514 and the G2 input data 
536 are interchanged, and 2) the red input data 516 and the Gl input data 538 are 
interchanged. In this configuration, the second-recovered color data on line 511 are 
G2 A data, the second-recovered color data on line 533 are R A data, the second-recovered 
color data on line 549 are B A data, and the second-recovered color data on line 551 are 
Gl A data. 

Another version of the second color recovery module 122 operates as 
illustrated in Fig. 6. This version illustrates the abbreviated procedure for the second 
color recovery in a 4-color system. The same figure applies to both odd and even 
numbered columns. Interpolator 602 receives second color R data 618 and interpolates 
this data in the second dimension to provide second-interpolated red data. 

Similarly, interpolator 616 receives second color G2 data 620 and 
interpolates this data in the second dimension to provide second-interpolated G2 data. 

Switch 604 and switch 606 select different colors to provide as inputs 
to combiner 610. In this situation, combiner 610 forms the difference between the R 
and G2 data corresponding to each individual image sensor element according to: 

(35) R-G2 1 , or 

(36) R'-G2. 
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The difference is provided via difference channel 61 1 to filter 612. 

Filter 612 can be a non-linear filter that is provided to remove narrow 
peaks and valleys within the color difference signal 611. Combiner 614 combines the 
second color data selected by switch 606 with the filtered difference data 613. In this 
embodiment combiner 614 provides: 

(37) R-F(R-G2), or 

(38) G2+F(R-G2) 

where F() represents the operation of filter 612. As can be seen from equations (37) and 
(38), the second-recovered G2 data G2 A are provided by equation (37) and the second- 
recovered red data R A are provided by equation (20). These recovered color data are 
provided via line 617. 

Interpolator 622 receives second color B data 638 color data and 
interpolates this data in the second dimension to provide second-interpolated B data. 

Similarly, interpolator 636 receives second color Gl data 640 and 
interpolates this data in the second dimension to provide second-interpolated Gl data. 

Switch 624 and switch 626 select different colors to provide as inputs 
to combiner 630. Combiner 630 forms the difference between the B and Gl data 
corresponding to each individual image sensor element according to: 

(39) Gl-B\or 

(40) Gl'-B. 

The difference is provided via difference channel 63 1 to filter 632. 

In this case, filter 632 can be a non-linear filter that is provided to 
remove narrow peaks and valleys within the color difference signal 63 1 and combiner 
634 combines the second color data selected by switch 626 with the filtered difference 
data 633. In one embodiment combiner 634 provides: 

(41) Gl-F(Gl-B), or 

(42) B+F(G1-B), 

where F() represents the operation of filter 632. As can be seen from equations (41) and 
(42), the second-recovered B data B A are provided by equation (41) and the second 
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recovered Gl data G1 A are provided by equation (42). These recovered color data are 
provided via line 637. 

Examples 

5 

Example 1 
(Bayer Pattern With Three Colors) 

The present invention is applicable to two-dimensional array of image 
sensing elements, wherein the elements are arranged in a so-called "Bayer Pattern." 
30 The disposition of colors in a preferred variant of such pattern is shown in the following 

5 diagram: 
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1 5 In step 1 , each row is digitally processed to recover any colors that exist 

in that row, particularly by applying the color recovery methodologies described in U.S. 
Patent No. 4,774,565 (W.T. Freeman) for a one-dimensional array. For odd numbered 
rows, this means that the red and green data are interpolated, the color difference R-G 
is formed, and the color difference is filtered along the rows with a 1-D median filter 

20 to produce (R-G) m f_ row . Then, improved estimates are found for the missing colors 

according to the prescription: 
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At R pixels: G A = R-(R-G) mf 



row 



At G pixels: R A = G + (R-G) mf 



The same basic steps are employed at even numbered rows, but using the 
colors G and B. After treatment of both the odd- and even-numbered rows, this results 
in the following information: 
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.etc. 



The odd numbered rows are now missing only an estimate for the blue 



intensity at each element, and even numbered rows are missing only an estimate of the 
red intensities. Values for green have already been either measured or estimated at 
every element, so the recovery of green intensities is complete at this point. 

In Step 2, the still-missing colors are recovered by applying the color 
recovery methodology mentioned in Step 1, but in a column-wise fashion. 

Recovery begins with a one-dimensional interpolation of the R and B 
values in the column direction. 

For the odd-numbered columns, two color differences 5 (R-B) and (R-G), 
are formed. These are filtered in the column direction with a 1-D median filter, and the 
missing colors are recovered as follows: 



At RG A elements: 



B A = R-(R-B) mLcoI 



At B A G elements: 



R A = G + (R-G) mf _coi 
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Even-numbered columns are treated similarly, except the required color 
differences are now (R-B) an (B-G). The missing colors are recovered according to: 
At R A G elements: B A = G + (B-G) m f_ CO i 
At G A B elements: R A = B + (R-B) mf _ C0 ] 

At this point, all colors have been recovered and the process is complete. 
Note that in all of the recovery operations, the recovered color at each element is 
5 computed as the combination of the one explicitly measured color with a filtered color 

difference. 

Example 2 

Z (Bayer-like Pattern with Four Colors) 

8 The method described in Example 1 is not restricted to a specific choice 

40 of colors, or to a specific disposition of colors. The following example illustrates this 

1 by application to a Bayer-like pattern in which there are two different shades of green 

(Gl and G2). This could happen by design, but in fact it is also likely to happen in 
3 practice even when a sensor is designed with only a single intended shade of green. 

7 An inspection of the Bayer pattern shows that there are two types of G 

j45 sites. The first G site (Gl) has R sensors located above and below, and B sensors 

a located to the right and left. The second G site (G2) has the B sensors located above 

and below, and the R sensors located to the right and left. Unless manufacturing 
methods are ideal, these two sites are likely to exhibit slightly different responses, and 
such differences may result in color fidelity problems. Even in the absence of 
20 manufacturing differences, the difference in structure of the sensing elements in the 

horizontal and vertical directions may lead to response differences due to electrical or 
optical "cross-talk" between pixels. Therefore, instead of treating these sites as 
equivalent, superior results is obtained by treating them as effectively different shades 
of green. The color layout is then: 



R Gl R Gl R Gl R Gl R Gl 
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G2 B G2 B G2 

R Gl R Gl R 

G2 B G2 B G2 

.etc. 



B G2 B G2 B 
Gl R Gl R Gl 
B G2 B G2 B 



In Step 1 of the recovery, each row is digitally processed to recover any 
colors that exist in that row, particularly by applying the methodologies described in 
U.S. Pat. No. 4,774,565 (W.T. Freeman) for a one-dimensional array. The result is: 



R,G1 A R A ,G1 

G2,B A G2 A ,B 

R,G1 A R A ,G1 

G2,B A G2 A ,B 

. .etc. 



R,G1 A R A ,G1 

G2,B A G2 A ,B 

R,G1 A R A ,G1 

G2,B A G2 A ,B 



R,G1 A R A ,G1 

G2,B A G2 A ,B 

R,G1 A R A ,G1 

G2,B A G2 A ,B 



In step 2, there are now two missing values per element. In the odd- 
numbered columns, for example, the R,G1 A sites are missing values for G2 and B, 
while the G2,B A sites are missing values for R and Gl . To remedy this, all four color- 
values are interpolated in the column-wise direction, and the required color differences 
(R-G2), (R-B), and (G1-G2) are formed. Then the missing values are estimated as: 



At R,G1 A sites: G2 A = R - (R-G2) mf co i 

B A = R - (R-B) mf col 

At G2,B A sites: G1 A = G2 + (Gl-G2) mf< 

R A = G2 + (R-G2) mfco , 



A similar pattern is followed for even-numbered columns, but using 
color differences (G1-G2), (B-Gl) and (R-B). 
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The procedure just described produces results with fidelity equivalent 
or superior to that of the two-dimensional form of color recovery described in U.S. Pat. 
No. 4,774,565 (W.T. Freeman), and the computations are considerably faster and 
simpler to implement in special-purpose hardware. 

5 

Example 3 

Examples 1 and 2 above refer to rows and columns of pixels, indicative 
of a rectangular array of sensing elements. It will be appreciated, however, that it is not 
a requirement of the present invention that the rows and columns be disposed at right 
40 angles to one another. For example, the present invention can be performed on a 

=f hexagonal array of pixels as shown in FIG. 8, using "rows" and "columns" 10 and 12 

J that are separated in angle by 60 degrees. Thus, the step of recovering missing color 

j 

5 information along a first dimension according to the invention is conducted, for 

Z example, along dimension of 1 0 of the hexagonal array, and the step of recovering 

1 5 missing color information along a second dimension is conducted along dimension 12. 

z The conduct of both is substantively equivalent to those identified in Example 1 . 

is Likewise, the results should also be substantively equivalent. 

While the invention has been described with reference to particular 
20 embodiments, it will be understood that the present invention is not limited to the 

particular constructions and methods herein disclosed and/or shown in the drawings, but 
also comprises any modifications or equivalents within the scope of the claims. 
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